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ABSTRACT 

Some of radiopulsars have anomalous braking index values n = VtVl/tl 2 ~ ±(10 3 ±10 4 ). 
It is shown that such anomalous values may be related to nondipolar magnetic field. 
Precession of a neutron star leads to rotation (in the reference frame of the star) of 
its angular velocity f2 around the direction of neutron star magnetic dipole moment 
m with an angular velocity Sl p . This process may cause the altering of electric current 
flowing through the inner gap and consequently the current losses on the time scale 
of precession period T p — 2tt/Q p . It occurs because the electric current in the inner 

gap is determined by Goldreich- Julian charge density pcj — that depends on 

the angle between direction of small scale magnetic field and the angular velocity fl. 
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1 INTRODUCTION 

Radio pulsars have been discovered more than 40 years ago 
I He wish! (|l9 68l and at present many thousand papers are 



devoted to these objects. Despite of large progress in the 
understanding of processes in pulsar magnetospheres still 
many important questions are unclear. One of them is the 
value of pulsar braking indices. If a pulsar were just a sim- 
ple magnet with dipolar magnetic momentum fn rotating 
with angular velocity fl, then the braking index n — Q.Q./tl z 
would be equal to 3. The taking into account of evolu tion of 
angle x between vector s m and Q yields n — 3 + 2cotx Davis 
l|l970l ); iBeskinl |2006al h Observations of some pulsars are 
well within theoretical predictions. For example, the brak- 
ing index o f pulsar Crab is equal to 2.5 and pulsar Vela 
has ra ~ 1. 4 lBeskinl (|2006al ). So it seems that more sophisti- 
cated magnetospheric mo dels will be able t o remove this dis- 
crepancy, see for exampl e Mclatos (1997); Timokhin (2005, 
2006); IContopoulosl (2006, 200^); lTimolAml |2007al lbl). How- 
ever, many isolated pulsars have very large positive or, some- 
times, negative braking indices upto |n| ~ 10 4 . Sometimes 
such larg e value s may be due to unobser ved glit c hes o r tim- 
ing noise lAlpail (|2006l V For example, in lAlparl |200d ) it is 
shown that all negative braking indices may be associated 
with unresolved glitches which has occurred between inter- 
vals of observations of the pulsars. But in some cases re- 
fined measurements of braking indices show that at least 



some pulsars have large a nd somet i mes n egative braking 
indices n ~ ±(10 — 10 ) IJohnstonl (|l999h . For example, 
the pulsars B0656+14 and B1915+13 have braking indices 
n ~ 14.1 ± 1.4 and n ~ 36.08 ± 0.48, correspondingly, 
while pulsar B2000+32 and B1719-37 have very large neg- 
ative braking in dices n ~ —226 ± 4.5 and n « —183 ± 10, 
correspondingly IJohnstonl |l999h . There are some explana- 
tions of such braking index values. For example, the val- 
ues like n ~ ±10 may be related to nonstandard mech- 
anisms of pulsar brakin g, like neutr on star slowing down 
due to neutrino emission iPenel (|l982p or due to interaction 
with c ircu mpulsar disk iMenorJ (|200ll ); iMalovl (|2004h ; IChenl 
l|2006l ). see lMalovl | 200 ll ) for review of possible mechanisms 
and its comparison with observations. The large braking in- 
dices may be explained by the rapid changing of magnetic 
fiel d. Such changes m a y be caused by Hall-d r ift ins tabili- 
ties iRheinhardl |2002r i; rGeppertJ (|2002l ); iPonsI l|2007h . Also 
the positive braking indices of old pulsars may be explained 
by the relaxation of angular velocit ies of neutro n star crust 
and superfluid between two glitches [Alparl (|2006l ). The large 
braki ng indic e s of so me pulsars may be related to Tkachenko 
waves IPopovl d2008l). In this paper we present a some m odel 
base d on work s TBeskinl (T2006bl ) and l ContopoulosI (120071 ). see 
also iBirvukwl (|2007l ); lUramar^Ood ); IPonsI (|2007r i. where it 
has been shown that large braking indices may be explained 
by the existence of some internal cyclic process in a pulsar. 
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2 THE TORQUE 

There are ma ny mechanisms , that may be responsible for 
pulsar braking [Beskml (|2006d ); lMalovl (|200lh . Let us shortly 
describe two of them. T he first of the m was proposed be- 
fore pulsar discovery by iPacinil |l967:). According to this 
mechanism the rotation energy E and the rotation mo- 
mentum M of a neutron star are carried away by mag- 
netic dipole radiation, that, of course, leads to a torque 
acting on the star. Because of the presence of currents and 
charges in pulsar magnetosphere, calculation of the strength 
of this torque is a highly complic ated and, a t pre s ent, un - 
completed t ask, s ee for e xample iTimokhinl (|2006l , l2007al ): 
ISpitkovskvl (|200d l; lBeski"nl |2006a). However, the torque can 
be estimated within a pure vacuum model of pulsar mag- 
netosp here, where ma gnetospheric currents and charges are 
absent IPeutschl l| 19551 ). In this case, the rotati on torque Kdi v 
acting on a ne utron star, can be written as iDavia l|l97d ); 
iMelatos! |2000l ): 



where jjv is the density of the electric current, that flows 
in the northern pulsar tube, and js is the density of elec- 
tric current that flows inside the southern pulsar tube, 
3gj C0S X — cos\ is the Goldreich- Julian current, Bo = 
2m/a 3 is the magnetic field strength at the magnetic pole, 
SN(rj) and Ss(ri) are the areas of the northern and south- 
ern pulsar tubes, So(rj) = iva 2 (Q,a/c)rj > , r\ = r/a, r is the 

distance from the center of t he s t ar. 

Following IJonesI (|l976l ); Ixiil (|200lD ; Iwiil i|2003l ), we as- 
sume that the star can slow down by both mechanisms si- 
multaneously and that the resulting torque K can be de- 
scribed just as a sum of partial torques Kdip and K cur : 



K = K dip + K c 



(6) 



Thus, the equation of rotation momentum loss can be writ- 
ten as 



dM 
~dt 



(7) 



Kdip = K (e m cos x — en + Rdi P [en x e„ 



(1) 



where H = S7e*n is the angular velocity of the star, m = me m 
is the dipolar magnetic momentum of the star, en and e m 
are two unit vectors, directed, correspondingly, along and 
m, x is the angle between en and e m . In the case of a small 
rotation speed Q,a/c <C 1, where a ~ 10 km is t he ra dius of 
the star, Ko and Rdi V can be written as IPavisI l|l970l ) 



2 m 2 ft 3 
K o = tt — ; — • 



Rdi, 



2 s Vfia 



cosx- 



(2) 



(3) 



In lDavisI l|l97Cl ): iGoldreicbl (jl970l ) the coefficient £, is taken 
to be equa l to 1. We guess that the value £ = 3/5 adopted 
in IMelatos! (2000) may be closer to reality, because it corre- 
sponds to absence of electric current sheet on the neutron 
star surface. Thus, we will further consider £ = 3/5. Expres- 
sion HI is valid at any values of Qa/c, although coefficients 
Ko and Rdi P slightly distingu ish from their nonrelativistic 
values, defined by © and (p)) . iMelato s! |2000l ). 

The second mechanism (current losses) is due to the 
electri c current j, that flow s along open f ield lines, see for ex- 
ample Heskin| (|l984 l2006al . 1 19931 , l2006cl ) . This current flows 
from light cylinder, crosses pulsar polar cap diode, intersects 
the star surface, and goes into deep crustal layers. After this 
it begins to move to surface and then transforms to the back- 
ward current flowing along the separatrix between open and 
closed field lines. It is worth to note that when the electric 
current travels inside the crust it sometimes flows across the 
magnetic field, that leads to torque K cur acting on t he star 
polar cap. The strength of this torque is calculated in I Jones! 
(j 19761) as: 



K cur = —Kq ■ a ■ e m cos x, 



(4) 



where parameter a(£l, y . </>) characterizes the electric cur- 
rent value Bcskin (2006c.a|). In the case of a circula r pulsar 
tube cr oss section t he pa rameter can be estimated as lBeskinl 
l|2006al lc3); [T"svganl (|2009l ) 



4 cosx 



3n ( S N (v) 
3°gj V 



3 SPHERICAL SYMMETRY CASE 

At first case we assume that a neutron star is an absolutely 
rigid sphere. Particularly, we neglect any star deformations 
and any viscosity and dissipation inside the star. Conse- 
quently, we can suppose that the star rotation momentum 
M = iCl, where / is the momentum of inertia of the star. 
Under this assumption equation (0 can be rewritten as 

I ~j]~ = K = Ko-( e m (l-a(f2,Y,0))cosx 

— en + Rdi P [en x e m ]). (8) 

Now let us to introduce three orthogonal unit vectors e x , 
e y and e z which rotate together with the star. As vector e 2 
we use e m , direction of vectors e x and e y may be arbitrary. 
Consequently, we have 

f =[«xe,],f =[iJx*],f§ =[Oxe H ] (9) 
and at any time t the following relations are valid: 



(e z ,e x ) = (e z , e y ) = (e x ,e y ) = and e 2 z = e\ 



1.(10) 



Hence, at any time t we can treat these vectors as an or- 
thogonal space basis. So it is possible to write, for example 



Q — Q (sin x cos <t>e x + sin x sin cf>e y + cos xe-n 



(11) 



Using the last expression, the equation (8} may be rewritten 
in the form 



(sin 2 x + a(n, X, <P) ■ cos 2 x) 
dx 1 

— = (1 - a(f2,x, </>))• sin x cosx 

at t 



dt t 
where 
~ 3 Ic 3 



Rdip — „ , 



2 m 2 !! 2 
2tv 



1.5 • 10 8 year 



(la) B?,' 



R 



~ _ 47r 1 fVla\ ~ 



■dip 



3 4 cos X V c / 

1.2 ■ 10 5 year ( -^—) ( £■) * 
V C cos x / Vis/ Bi 



2 i 
12 



(12) 
(13) 
(14) 

(15) 
(16) 
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where B12 = Bo/10 12 G. It is worth to note that equations 
(|12H14[) describe neutron star rotation in an inertial reference 
frame of "rigid stars" . But all vectors are represented by the 
basis {e x , &y, e*z) altering in time. 

It is easy to see that equation (|12|) describes the losses 
of neutron star rotation energy and equation (|13[) describes 
the evolution of the inclination angle X - I n the case of a = 1 
the equation (|12|l may be rewritten as 



dP _ P 
dt ~ f' 



(17) 



So in this case parameter r is equal to two characteristic 
pulsar ages r = P/{2P). For any other values of parameter 
a there is no so simple interpretation of the time r. But it is 
possible to treat the parameter r just as some characteristic 
time over which large changes of the angular velocity Q. and 
the inclination angle x occur. The equation (|14|) describes 
the precession of the neutron star and T v may be interpreted 
as precession period. It is worth to note that the period T p is 
3-4 orders of magnitude smaller than the characteristic time 
t. So it seems that the neglecting of the altering of values 
f2 and \ over times compared with precession period T p 
may be a good approximation. Consequently, if we assume 
a = a{Q, x, 4 1 ) h will be possible to write 



da 
~d~t 



da d(f> 
d<f>~dt' 



(18) 



In this case, in equation (|12p only parameter a is able to 
change significantly over precession period and hence 



d 2 n 

dt 2 



da 2 

T' dT C0S x 



(19) 



Thus, the b raking index n — flQ/tl 2 can be estimated as 
iGoodi (|l985h 



an 



d(p da 
~dt ' ~b~4> 



cos 2 x 



(sin x + a cos X) 



(20) 



If we take into account the equation (|14|l . expression (f20 
may be rewritten as 



n « R, 



8a 



dip 



cos x 



30 (s 



sin x + a cos X) 



(21) 



If we assume that a ~ 1, cos x ~ sin X ~ 1> then the braking 
index n may be estimated as 



It J^dip ' r\ 1 • 

0<p 



(22) 



If we take into account the equation © then the braking 
index estimation (1211) takes the form 



2 s \Qa 



da 
~d~4> 



cos 2 x 



(sin 2 x + a cos 2 x) 



VfiJ d<h \P) 



da 



(23) 



It is necessary to note that within the current paper we 
assume that the magnetic dipole moment m does not change 
significantly over the precession period T p . Particularly, the 
expression (|23p is valid only if the time corresponding to 
large changes of the magnetic dipole momentum rh is large 
enough compared with the period T v . 



Now following iGoldreichl (|l97Ch : iTsyganl (|2009l ) we de- 
fine the value of function f{Q, x, <t>) averaged over precession 
time as 

1 f 2n 

<f>{Q, x ) = — J /(Q, X , <!>)d<t>. (24) 
Then, averaged equations (|12|) and (|13l) may be written as 

(25) 
(26) 



d^l £~2 / • 2 2 \ 

< — (an x+ < « > cos xj 

< — >= — — (1— < a >) sinxcosx- 



It is easy to see that case < a >= 1 corresponds to "station- 
ary" or, better to say, equilibrium value of the inc l inatio n 
angle, when < ^ >= o lGoldreicbl l|l97d ): lGure : dcn1 (|2007h : 
llstominl (|2007m . In this case the average speed of angular 
velocity decrease also does not depend on the inclination 
angle: 



n 



(27) 



dQ. 

The las t equation may c orrespond to the results of iBeskinl 
l|2006bD ; lBirvukovl ([2007) where it is shown that an average 
braking index of many pulsars is close to n ~ 5. It is worth to 
note that if the parameter < a > decreases with inc r easing 
angle %, this equilibrium state is stable, see ITsyganl (120091 ) 
for details. 



4 NONDIPOLAR MAGNETIC FIELD 

It is widely accepted that besides large-scale dipolar mag- 
netic field, small scale nondipolar magnetic field may exist 
near neutron star surface. This field has a spatial scale about 
i ~ (0.3 — 3)km and rapidly falls to zero when the distance 
from the neutron star surface increases to become fully neg- 
ligible at the light cylinder. This allows to suppose that this 
component does not exert any direct influence on magnetic 
dipole braking, that appears to occur near the light cylinder. 
It seems that the same reasons are applicable to the direct 
influence of the small s c ale co mponent on current losses, 
see for example ITsyganl (|2009| ). However, the existence of 
the small scale component may lead to drastic changes of 
the ele ctric cur r ent th at flows through the inner gaps of the 
pulsar IShibatal (|l99ll ). Thus, the indirect influence of the 
nondipolar field on current losses and torque K cur is sup- 
posed to play a decisive role in pulsar braking. 

The influence of small scale magnetic fiel d upon electri c 
cur r ents was considere d in many papers, c.f. lAronsI <|l9TSh : 



[GU l|200ll . 2003, 200630) ■ 

In order to estimate its influence 

on curr e nt los ses we wi ll use a simpl e model, propos ed in 
IPalshinl (| 19981 ). see also iKantorl (|2003l ): ITsyganl (|2009l ). Ac- 
cording to this model the inner gaps fully occupy the pulsar 
tubes cross sections and are situated close to the neutron 
star surface. It is assumed that magnetic field strength on 
the star surface is not large enough to prevent free emis- 
sion of electrons, so when pulsar diode is placed onto star 
surface it will be operate in electron charge limited steady 
flow regime. In the neighbourhood of the inner gap small 
scale magnetic field will be modeled by the field of a small 
magnetic dipole fni, which is embedded into neutron star 
crust under magnetic (dipolar) pole at the depth A • a (see 
fig. [T] and fig. [2]). We assume that A w 1/10 and that mi 
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is perpendicular to the main dipole momentum m. So the 
vector mi may be written in the form 



mi = 2mv/S. (e x cos 7 + e y sin 7) . 



(28) 



where v — B1/B0 is ratio of small scale magnetic field 
strength Bi = mi/(A ■ a) 3 to the strength Bo = 2m/ a? 
of the large scale dipolar magnetic field at the magnetic 
(dipolar) pole. We also assume that the changes of values v, 
A and 7, as well as, the value of magnetic momentum m, 
are negligible, at least, over a period T p of neut ron s tar pre - 
cession. A similar model has been proposed in Icil (|2002al l . 
where a small dipole rn 1 is placed not strictly under the 
magnetic pole and may have arbitrary direction. The case 
of pure axisym met r ical sm all scale magne ti c field has bee n 
investigated in[Gi3 (|2002al lbh; [Assed (|2002l ); iTsyganl (|2000D . 
It is worth to note that the field of dipole rh\ is able to 
describe small scale magnetic field only in vicinity of the 
pulsar diode. We does not any intention to suppose that 
this dipole is able to describe small scale magnetic compo- 
nent over whole neutron star surface. We rather suppose 
that the small scale field is close to a sum of 10 2 — 10 3 such 
small dipoles and that we use only one of them, which is the 
closest to pulsar diode. 

In the case of a thin pulsar tube R± <C A ■ a, where Rt is 
its radius, and at small altitudes 77 = r/a <C k^J c/Sla ~ 10 2 
the Goldreich- Julian density (inside the pulsar tube) may 
be written as 



pGJ — 



SIB 

2m 



(29) 



At v < 1/3 the function f(j}) may be approximated as 
IPalshirJlj 19981 ) 



cosx 



+A 



1 + g-l! ) sinxcos^- 7) 



(30) 



where the coefficie nt k ~ 0.1 5 descr ibes general relativis- 
tic frame dragging iMuslimovl (119921 ) . B - magnetic field 
strength, A = v (Ar?) 3 / (77 — 1 + A) 3 . In the case of k = 
the function f(rj) is equal to the cosine of the angle between 
the direction of magnetic field B and angular velocity SI. 

In the case of a thin long inner gap, when the radius Rt 
of pulsar tube is small compared with the gap height z c a, 
see fig. [TJ the density j of elect ric current flow ing through 
the inner gap may be written as I Tsygan] (|2009l ) 



3 « JajfiVo), 



(31) 



where 770 is the altitu de wher e the inner gap begins and 
3%j = §f . Following rTsvganl l|2009l ) we assume that the 
electric potential $ monotonically increases inside the pul- 
sar diode and the diode is placed as close to the surface 
as possible. Hence using expression (|30|) one may write at 
cos(0 — 7) < 



j(y,i) « Jgj - max (0, /(l)) 
and in the case of cos(<^> — 7) > 
iiy, 7) w 3gj ■ min (/(??), 1 < rj < +00) , 
see lTsyganl (|2009r i for details. 



(32) 



(33) 



northern inner gap may be described by a small dipole with 
parameters v — vn and 7 = 7jv, and the small scale mag- 
netic field in a neighbourhood of southern inner gap (and, 
of course, inside it) may be described by a small dipole with 
v = us and 7 = 75. Then, the parameter a may be written 
as 

( 3\vn,1n) ( S N (y) \ KysTfs) ( S s (y) \ \ 



4 cosx 

where Sn and Ss is cross section areas of northern and 
southern pulsar tube. 

Let us firstly consider the most simple case when i>n = 
vs = v and Sn — Ss — So- The values of braking index n 
for angle x = 30° are shown on fig. [3] and [3] They have been 
calculated by the substitution of expressions (|30[1 and ((34} 
into equation (|21[) . It is easy to see that the values v ~ 0.1 
mostly correspond ton~ 10 — 10 2 and for v ~ 1 the braking 
index may be as large as ~ 10 3 . 

On fig. [7] and |8] it is shown the braking indices when 
angle x is not arbitrary but correspond to equilibrium value 
Xeq, which is defined as < a > (x = Xeq) — !• It is easy 
to see that in this case at v = 0.1 the braking index n ~ 
±(10 — 20) and at v — 0.7 the braking index may achieve 
values as large as n ~ 1.5 • 10 3 . 

In the case of weak small scale magnetic field v <C 1 
and viaxix < 1 the electric current through the inner gap 
may be approximated as 



Jgj cos X ( 1 - « - v tan x > 



•7)0(008(^-7))), (35) 



where O(x) is the Heaviside function, <d(x) — 1 at x ^ 
and Q{x) = at x < 0. Consequently, parameter a is equal 
to 

3 1 



(36) 



4 cos x 3gj 

The value of parameter a averaged over precession period 
T p may be written as 

3 
2 



< Q > 



VN + VS 



2tt 



tan x 



(37) 



Thus, the equilibrium angle x is equal to 
' 2tt 1 

v l/jv + vs 3 



(38) 

At this angle the expression (1361) may be rewritten as 

- 3k) x 

- 7jv)6(cos(<?!> - 7at)) 

>-7s)e(cos(0-7s)))) (39) 
vs and 7jv = 7s = parameter a 



+ 



vn 
vn + vs 
vs 



VN + VS 

and in the case of vn 
may be estimated as 



3 / 1 

— ( 1 — — cos (f> ■ G(cos c 



(40) 



This expression shows that at the equilibrium value of angle 
X parameter a is changing from w | to w | or, in other 
words, is changing in two times over the precession period. 
Consequently, the braking index n may be a crudely esti- 
mated as n ~ 2ttt/T p 



10 3 



Suppose that the small scale magnetic field near the 



Now let us discuss a more sophisticated case when the 



The braking indices of radiopulsars 5 




Figure 1. The position of small dipole mi is shown in the case 
of 7 = and = 0. The lower gray square is neutron star and 
upper dark area is pulsar diode. The pulsar diode is situated on 
star surface t)q = 1, Rt is the radius of pulsar tube and z c a is the 
height of the diode. 




Figure 2. The orientation of vectors Q, m and mi in the case 
of 7 = 0. 



pulsar tube cross section St depends on angle y. We will use 
the dependence that was derived in iBiggsj ([l990): 



St(v) « S (v)g 2 (x) 

where 



g(x) = 



(1 _ f) 3xl/4 



(41) 



(42) 



and fj, = 1 + § cos 2 x 



v i cos xl \/2 + \ cos 2 x- When % in- 
creases from to -| the value /j, changes from to 1/3. At 
X — 0° the coefficient g(x) is equal to 1 and at x — 90° it 

is equal to g (jj = (57)* ~ 0.620, so the area of pulsar 
tube cross section decreases with increasing angle X- Some 
other dependencies of puls ar tub e area s on the an g le y i s 
calcula t ed, for example, in iDvksl (|2004T ); iMuslimovl {2009); 
iBeskinl (|2006ah . Such a depende nce may be eit her decreas- 
ing or increasing, for example, in lBeskinl l|2006al ) it is shown 
that the area of a pulsar tube increases with x from 1.59So 
at x = 0° upto 1.96S at % = 90°. 

In the case of vn = vs = v an d Sn — Ss = Sog 2 (x) 
the resulting braking indices n for angle x = 30° are shown 
on fig. [S] and [5] Braking indices at equilibrium values of the 
angle x are shown on fig. [9] and 1101 



5 AN AXISYMMETRIC CASE 

In this section we will suppose that a neutron star is an ab- 
solutely rigid axisymmetric body and that its symmetry axis 
coincides with dipolar magnetic momentum rfi. The rotation 
momentum of the star may be written as 



M 



E 

0=x 



Ip-ep- (e> ■ n) = I ± Q - AJe m (e m ■ Q) , (43) 



where I x = I v = I± 



and I z = In = I a 



— AI are momenta 



of inertia of the star. Hence, the equation (J7) of momentum 
loss takes the form: 



(44) 




Figure 3. The dependence of braking index n on precession 
phase (f>. The changing of phase <f> from — n to ir corresponds to 
one precession period and takes T p seconds. It is assumed that 
S N = Ss = So, vn = vs = v, 7iv = 7S = 0, x = 30°, P = Is 
and the neutron star is spherical. 



where K e f / is an effective rotation torque acting on the star 

K e ff = K (e m (1 - a(Q, x, </>)) cos X - eh 

+Reff[en x e m ]) . (45) 

Here we neglect the small term AIe m ^e m • = 
m (e m • K) and introduce the coefficient 



(A/// Z )e„ 



,Am 2 



2 \naj- 



cosx 



AIc 2 a 



cosx 



3- 10 



(£)(! 



aj 33 

13 12 



cosx, 



(46) 



where A7 33 = A7/10 33 g cm 2 . 

It is easy to see that the only difference between equa- 
tions ([8j and (|44|l is the replacing of momentum of inertia 
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Figure 4. The same as fig. [3] but 7jy = and 73 = 



Figure 7. The dependence of braking index n on precession 
phase (f>. The changing of phase <f> from — tt to 7r corresponds to 
one precession period and takes T p seconds. It is assumed that 
Sn = S s = S , v N = vs =v,1N =1S = 0, x = Xeq, P = Is and 
the neutron star is spherical. Equilibrium angles Xeq are equal to 
X eq ~ 81° at v = 0.1, Xeq ~ 65° at v = 0.3, Xe 9 « 49° at v = 0.5 
and Yen ~ 33° at v = 0.7. 
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Figure 5. The same as fig. [3] but Sjv = &S = •Sb<? 2 (x)- 
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Figure 6. The same as fig. but 7jy = and 73 = 



Figure 8. The same as fig. but 7jv 
of equilibrium angles Xecj is the same. 



and 7s = . Values 



/ by its component I± and the replacing of the coefficient 
Rdip by the coefficient R e ff- Consequently, upto these two 
exchanges all the formulas of previous section remain appli- 
cable to the axisymmetrical neutron star. Particularly, the 
braking index n may be estimated as 



n « R, 



da cos 2 X R da 

!// '^'( si n 2 X + aco S 2 X ) 2 ~ Heff 'd4>- ' 



Again, supposing that a ~ 1, cosx ~ sm X 
crudely estimate the braking index as 



1, one can 



Reff sin 

,4 (P 
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Figure 9. The same as fig. [7] but Sjv = Ss = Sog 2 {x)- Equilib- 
rium angles \eq are equal to \eg ~ 46° at v = 0.1, \eq ~ 39° at 
v = 0.3, Xeq ~ 30° at v = 0.5 and \eq ~ 21° at u = 0.7. 



Figure 11. The dependence of braking index n on precession 
phase <j>. The changing of phase <f> from —n to 7r corresponds to 
one precession period and takes T p seconds. It is assumed that 
Sn = S s = S , vn = v s = v, 7jv = 7S = 0, x = 30°, P = Is, 
B = 10 12 G and AI = -10~ 13 I. 
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Figure 10. The same as fig. [9] but 7jv 
of equilibrium angles Xeq is the same. 
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Figure 12. The same as fig. El but AI = -2 • 10~ 13 /. 



250year 



(48) 



where r = P/ (2P) is characteristic pulsar age and J ~ Jx is 
moment of inertia of a undeformed (spherical) neutron star, 
7 45 = //10 45 gcm 2 . 

Some examples of dependence of braking index n on the 
precession phase (j> are shown on fig. I1U13I This dependence 
may be found by calculating the braking index n BV h £r for a 
spherically symmetrical neutron star and then multiplying 
it by R e ff/Rdi P : 



eff 



Rdi 



i^spher • 



(49) 



The axisymmetrical deformations of the neutron stars 
may be caused by internal magnetic fiel d that res i des in - 
side neutron star crust or inside its core iGoldreichl (|l970h . 
The deformation caused by internal magnetic field may be 
estimated as 



A7 
I 



r>2 A 



(50) 



GM 2 ' 

where Bi„ is the strength of the internal magnetic field and 
M is the mass of the star. The coefficient £ depends on 
magnetic field profile. In the case of dipo l ar ma gnetic field 
£ may be estim ated as C = 25/8 iFerrarol (I1954T ). the same 
value is used in IGoldreichl (|l970l ). For othe r configuration s 
it may be, for example, as small as £ = 1/18 Haskell (2008). 
If we assume that B in — Bo = 2m/a 3 then the coefficient 
Reff can be estimated as 



R. 



eff 



' r g Ma 2 



cos x 



G£)H 



3- 10' 



(51) 



where r g = 2GM/c 2 is gravitational radius of the neutron 
star. Firstly it is easy to see that in the case of £ ~ 1/18 the 
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Figure 13. The same as fig. [TT] but AI ■ 
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Figure 14. Dependence of value / = j| n (^p) on the char- 
acteristic pulsar age r = ^ for various pulsars. The solid line 
corresponds to value / determined by l|48j l and (1 50 1) f = 25/8 and 
the dashed l ine co rresponds to ( = 1/18. Pulsar data is taken 
from lATNlj l|20ldV 



braking indices does change significantly. Secondly if £ > 1 
then the braking indices will be ~ 20 • £ times more than in 
the spherically-sy mmetrical case a nd, particularly, at value 
C = 25/8 used by iGoldreichl (|l970h the braking index could 
become as large as n ~ 10 4 - 10 s . 

In some cases the neutron star interiors may be super- 
conductive. If the neutron star matter is a superconductor 
of the second type, the internal magnetic field is able to form 
magnetic flux tubes. In this case the coefficient £ increase 
substantially £ ~ B f i/B in ~ 10 3 , where B n ~ 10 15 G is the 
strength of magnetic field inside magnetic flu x tubes and Bj n 
is average strength of internal magnetic field lAkgunll|l987f) . 
This leads to Bfi/Bi n ~ 10 3 time larger deformation of the 
star and, consequently, to very short precession period like 

P P ~ — (pi years . (52) 
cosx V Is/ 

Hence, the coefficient R e ff increases significantly and leads 
to the increasing of braking index, that can reach n ~ 10 7 — 
10 s . As the absolute majority of normal isolated radiopulsars 
do not have such large braking indices it may be possible to 
conclude that there is no superconductivity of the second 
type in neutron star interiors. 
The dependence of value 

'-Mt) 

on characteristic p ulsar age r is s hown on fig. 1141 The pulsar 
data is taken from lATNFI (|2010h - If the estimation ((48j) of 
braking index and the estimation (|51[1 of neutron star de- 
formation were valid then value / would be depend only on 
parameter (. The solid line shows the estimation of value 
/ at £ = 25/8 and dashed line corresponds to the case of 
£ = 1/18. Because of the estimation (|48p yields only the 
upper limit of possible values of braking indices the obser- 
vation data are in good agreement with the case of £ = 25/8. 
Although the majority of pulsars does not contradict to the 
case of = 1/18 too. 

On fig. [15] it is shown the dependence of value / on 
dipolar magn etic field Bo- The pulsar data is taken from 
ATNl] (|2010l ). Tho horizontal lines correspond to constant 
values of parameter Solid line corresponds to £ = 25/8 and 
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Figure 15. Dependence of value / = §"-(^r) on the dipo- 
lar magnetic field Bq for various pulsars. The solid line corre- 
sponds to value / determined by (I486 and (1501 ) £ = 25/8 and 
the dashed line corresponds to £ = 1/18. The two dot-dashed 
lines correspond to the case of ( = ^-{Bji/Bq) (upper line) and 
C = j %(Bfi/Bo) ( lower line), B/i = 10 15 G. Pulsar data is taken 
from lATNFI teOloT) . 

dashed line corresponds to £ = 1/18. The two dot-dashed 
lines correspond to the case of C, — (25/8) • (Bfi/Bo) (upper 
line) and ( = (1/18) • (B fl /B ) (lower line), B fl = 10 15 G. 
The last two lines correspond to the case of neutron star 
deformations caused by magnetic flux tubes and the average 
internal magnetic field Bi„ is equal to Bo- 

The formula (|48[) allows to estimate the lower limit of 
neutron star deformation AI/I. The corresponding values 
are shown on fig. 1161 It shows that braking indices of ma- 
jority of pulsars may be explained by the existence of defor- 
mation AI/I ~ 1 0~ 13 — 10~ n tha t agrees with estimation 
of AI received by IGoldreichl jl970l ). 

6 CONCLUSION 

In this paper we present a some explanation of large values 
of braking indices of pulsars. The proposed model is based 
on four main assumptions: 

(i) The neutron stars have a small scale magnetic field. 
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Figure 16. The estimation of a lower limit of neutron star defor- 
mation AI/I by braking indices for various pulsars. Value AI/I 
is estimated as K. %n ( — ) r m ,) . Pulsar data is taken from 
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precession periods T p larger than (10 2 — 10 4 )P may exist 
only when vortices of superfiuid can not been pinned any- 
where in star, are a ble to move freely a nd do not coexist with 
magnetic flux tubes [Lmkl (2006, 20q3)- The small number of 
isolated pulsars with observed precession force us to exclude 
the triaxial precession and to assume that neutron star is 
axisymmetrical and symmetry axis coincide with magnetic 
dipole moment rh. In this case in frame reference related to 
"rigid stars" vector m rotate with constant angular velocity 
uj. And because of uj ~ Q its trajectory coincides with the 
case of unprecessing star, see Appendix A. Also, in order to 
prevent the observation of such precession it is necessary to 
assume that pulsar tube structure does not precess too. It 
particularly means that pulsar tube crossection is circular 
or depends on only vectors Q, and m and does not depend 
from small scale magnetic field. Also, it means that distri- 
butions of energy of primary electrons and pair multiplicity 
over pulsar tube crossection are close to axisymmetrical. 



The strength of this field must be enough large to curve pul- 
sar tube but enough small to allow free emission of electrons 
from star surface. 

(ii) There are inner gaps in pulsar tubes and a electric 
current flowing across inner gaps depends on angle between 
small scale magnetic field and angular velocity Q of the star. 

(iii) Braking torque depends on current that flows across 
inner gaps. 

(iv) The neutron stars precess with periods like 10 3 — 10 4 
years. 

In this paper we also neglect the contribution of outer gaps. 
If a some part of electric current flows through outer gaps 
then this part is determined only by outer gap electrody- 
namics and, consequently, does not depend from small scale 
magnetic field. It decreases the variation of current losses 
over the precession period and, consequently, leads to the 
decreasing of braking index. In the case of young pulsars like 
Crab and Vela, which have small braking indices, we suppose 
that current losses is almost fully determined by currents 
flowing through outer gaps and the contribution of inner 
gaps current is negligible. In this paper it is also assumed 
that magnetic dipole braking exists and does not depends on 
electric current flo wing acros s inne r gaps . This as sumption 
is wid e ly used c.f. lEliseeval (|2006l ); lYuel (|2007l); iGurevichl 
l|2007f l: llstominl (|2007f) . but, as mentioned in lBeskinl (1200681 ). 
it must be treated with caution. It is shown that in the case 
of force free magnetosphere and absence of electric current 
flowing along pulsar tu be the orthogonal pulsar y = 7r/ 2 
does not slow down at all lBeskinl (|l983l . ll984l ); lMeste]| (|l999h . 
It gives the reason to suppose that magnetic dipole braking 
does no t exist o r, at le ast, must be depend on the electric 
current [Beskin (2006a). In such case presented model can 
provide the only qualitative explanation of the existence of 
large braking indices. And the quantitative estimations, of 
course, will strongly depend on relation between magnetic 
dipole braking and current losses torques. 

The presence of a long period precession is the weakest 
point of the model. At present the precession is discovered 
only at a few isolated neutron stars. And these pulsars hav e 
the precession periods like T p ~ 1 — 10 years iLmkl l|2007h . 
The presence of pined superfiuid in n eutron st a r crus t sub- 
stantially increases precession speed IShahaml l| 19771 ). The 



ACKNOWLEDGMENTS 

This work was supported by the Russian Foundation for 
Basic Research (project 10-02-00327) and the Program of 
the State Support for Leading Scientific Schools of the 
Russian Federation (grant NSh-3769.2010.2). We sincerely 
thank A.M. Krassilchtchikov for help with the preparation 
of the paper, A.I. Chugunov for help in the calculations 
and D.G. Yakovlev, A.I. Chugunov, M.E. Gusakov, Yu.A. 
Shibanov, D.A. Zyuzin, A.A. Danilenko, E.M. Kantor, V.D. 
Palshin and V.A. Urpin for help, comments, and usefull dis- 
cussions. 



REFERENCES 

T. Akgun and I. Wasserman, 2008, MNRAS, V. 383 p. 
1551. 

M.A. Alpar and A. Baykal, 2006, MNRAS, V. 372, p. 489. 

A. Alpar and H. Ogelman, 1987, Astronomy and Astro- 
physics, V. 185, p. 196. 

J. Arons and E.T. Scharlemann, 1979, ApJ., V. 231, p.854. 

E. Asseo and D. Khechinashvili, 2002, MNRAS, V. 334, p. 
743. 

ATNF Pulsar Catalogue, 

|http://www.atnf. csiro.au /research/pulsar /psrcat 

Beskin V.S. " Axisymmetric stationary flow in astro- 
physics" (Fizmatlitt, Moscow, 2006) in Russian 

G. Beskin, A.Biryukov and S.Karpov, arXiv:astro-ph/ 
0603375 

V.S. Beskin, A.V. Gurevich and Ya. N. Istomin, 1983, 
ZhETP V. 85 p. 401, in Russian 

V.S. Beskin, A.V. Gurevich and Ia.N. Istomin, 1984, As- 
trophysics and Space Science v. 102, p. 301. 

V.S. Beskin, A.V. Gurevich, Ya.N. Istomin, "Physics 
of pulsar magnetosphere", Cambridge University Press 
(1993) 

V.S. Beskin and E.E. Nokhrina, 2006, Astrophysics and 

Space Science v. 308, p. 569. 
J.D. Biggs, 1990, MNRAS V. 245, p. 514. 
A. Biryukov, G.G. Beskin, S. Karpov and L. Chmyreva, 

2007, Advances in Space Research, V. 40, p. 1498. 



10 D.P. Barsukov and A.I. Tsygan 



W.C. Chen and X.D. Li, 2006, Astronomy and Astro- 
physics, V. 450, p. 1L. 

I. Contopoulos, 2007, Astronomy and Astrophysics, V. 
466, p. 301. 

I. Contopoulos, 2007, Astronomy and Astrophysics, V. 
475, p. 639. 

I. Contopoulos and A. Spitkovsky, 2006, ApJ, V. 643, p. 
1139-1145. 

L. Davis and M. Goldstein, 1970, ApJ, v. 159, p. 81L. 
A.J. Deutsch, 1955, Annales D'Astrophysique v. 18, p.l. 
J. Dyks, A.K. Harding and B. Rudak, 2004, ApJ V. 606, 
p. 1125. 

S.A. Eliseeva, S.B. Popov and V.S. Beskin, arXiv : astro-ph 

/ 0611320 
V.C.A. Ferraro, 1954, ApJ, V. 119, p.407. 
U. Geppert and M. Rheinhardt, 2002, Astronomy and As- 
trophysics, V. 392, p. 1015. 
J. Gil, G.I. Melikidze and U. Geppert, 2003, Astronomy 

and Astrophysics, V. 407, p. 315. 
J.A. Gil, G.I. Melikidze and D. Mitra, 2002, Astronomy 

and Astrophysics, V. 388, p. 235. 
J.A. Gil, G.I. Melikidze and D. Mitra, 2002, Astronomy 

and Astrophysics, V. 388, p. 246. 
J. Gil, G. Melikidze and B. Zhang, 2006, Astronomy and 

Astrophysics, V. 457, p. 5L. 
J. Gil, G. Melikidze and B. Zhang, 2006, ApJ, V. 650, p. 

1048. 

J. Gil and D. Mitra, 2001, ApJ V. 550, p.383. 

P. Goldreich, 1970, ApJ, V. 160, 11L. 

M.L. Good and K.K. Ng, 1985, ApJ, V. 299, p. 706. 

A. V. Gurevich and Ya.N. Istomin, 2007, MNRAS, V. 377, 
p.1663. 

B. Haskell, L. Samuelsson, K. Glampedakis and N. Ander- 
sson, 2008, MNRAS, V. 385, p. 531. 

A. Hewish, S.J. Bell, J. D. H. Pilkington, P.F. Scott, R.A. 
Collins, 1968, Nature, V. 217, p. 709. 

Ya.N. Istomin and T.V. Shabanova, 2007, Astronomy Re- 
ports, V. 51, p. 119. 

P.B. Jones, 1976, ApJ, V. 209, p. 602. 

S. Johnston and D. Galloway, 1999, MNRAS, V. 306, p. 
50L. 

E. M. Kantor and A.I. Tsygan, 2003, Astronomy Reports, 
V. 47, p. 615. 

B. Link, 2003, "Radio Pulsars" ASP Conference Procs, V. 
302, p. 241. 

B. Link, 2006, Astronomy and Astrophysics, V. 458, p. 881. 
B. Link, 2007, Astrophysics and Space Science, V. 308, p. 
435. 

I.F. Malov, 2001, Astronomy Reports, V. 45, p. 389. 
I.F. Malov, 2004, Astronomy Reports, V. 48, p. 337. 
A. Melatos, 1997, MNRAS, V. 288, p. 1049. 
A. Melatos, 2000, MNRAS, v. 313, p. 217. 
K. Menou, R. Perna and L. Hernquist, 2001, ApJ, V. 554, 
p. 63L. 

L. Mestel, P. Panagi and S. Shibata, 1999, MNRAS, V. 
309, p. 388. 

A.G. Muslimov and A.K. Harding, 2009, ApJ, V. 692, p. 
140. 

A.G. Muslimov and A.I. Tsygan, 1992, MNRAS V. 255, 
p. 61. 

F. Pacini, 1967, Nature, V. 216, p. 567. 

V.D. Palshin and A.I. Tsygan, Preprint of Physico- 



Techinical Institute N° 1718, Saint-Petersburg (1998) in 
Russian 

Q.H. Peng, K.L. Huang and J.H. Huang, 1982, Astronomy 
and Astrophysics, V. 107, p. 258. 

J.A. Pons and U. Geppert, 2007, Astronomy and Astro- 
physics, V. 470, p. 303. 

S.B. Popov, 2008, Astrophysics and Space Science, V. 317, 
p. 175. 

M. Rheinhardt and U. Geppert, 2002, Physical Review Let- 
ters, V. 88, id. 101103. 

M. Ruderman and J. Gil, 2006, Astronomy and Astro- 
physics, V. 460, p. 31L. 

J. Shaham, 1977, ApJ, V. 214, p. 251. 

S. Shibata, 1991, ApJ, V. 378, p. 239. 

A. Spitkovsky, 2006, ApJ, V. 648, p. 51L. 

A.N. Timokhin, arXiv:astro-ph / 0507054 

A.N.Timokhin, 2006, MNRAS, V. 368, p. 1055. 

A.N. Timokhin, 2007, Astrophysics and Space Science, V. 
308, p. 575. 

A.N. Timokhin, 2007, MNRAS, V. 379, p. 605. 

A.I. Tsygan, "Pulsar Astronomy - 2000 and Beyond" ASP 
Conference Series, IAUS, V. 202, p. 473 (2000) 

A.I. Tsygan, P.I. Polyakova and D.P. Barsukov, 2009, As- 
tronomy Reports, V. 53, p. 1146. 

J.O. Urama, B. Link and J.M. Weisberg, 2006, MNRAS 
Letters, V. 370, p. 76L. 

F. Wu, R.X. Xu and J. Gil, 2003, Astronomy and Astro- 
physics, V. 409, p. 641. 

R.X. Xu and G.J. Qiao, 2001, ApJ, V. 561, p. 85L. 

Y.L. Yue, R.X. Xu and W.W. Zhu, 2007, Advances in Space 
Research, V. 40, p. 1491. 



APPENDIX A: THE PRECESSION OF 
AXISYMMETRICAL STAR 

Consider the neutron star rotation in an extreme case when 
it is possible to neglect the first and second terms in ex- 
pression (|45p which contain vectors e m and eh compared 
with the last term which contains [en x e m ] . Particularly, it 
means that any rotation energy losses are neglected and it is 
assumed that the angular velocity fi and inclination angle \ 
do not change with time. Consequently, this approximation 
is valid only over time scales that are small compared with 
the characteristic time r or pulsar age r, although this time 
scales may be comparable with or larger than the precession 
period T p . In this case the equation (|44[1 may be written as 

= ^Refffi x e m ]. (Al) 

With this equation it easy to obtain d^l/dt = and dx/dt = 
and, consequently, neither the characteristic time r nor the 
coefficient R e ff are changing. 

Let us introduce the vector e u = cos Pen — sin /3e m and 
calculate its time derivative 

_= _ (sm/3en+cos/?em) _ 

+ (_f« COS/ 0_fisin/3)[en x e m ] = 0. (A2) 
r 

Thus, in order to make vector e u constant in time, it is 
enough to choose the angle ft as: 
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Figure Al. The rotation of vector e a with angular velocity Q p . 

/3 = -arctan ( ^U-) . (A3) 
\ Sit / 

In this case the first of equations ([9]) may be rewritten as 

-rr = j^xem. (A4 

at cos p 

It means that e m and, consequently, m just rotate with 
the constant angular velocity Co = eLfl/ c os/3 a round the 
constant vector Co, see for example [LmM (|2003h . As Q = 
Q tan(/3), would also rotate with the angular ve- 
locity Co around the constant direction e u . It is worth to note 
that the angle /3 is very small: 

P = arctan + cos X 

and, consequently, Co is almost undistinguished from the star 
angular velocity f2. It seems that such a rotation substan- 
tially increases the difficulty of direct observations of pulsar 
precession. It is easy to see, that in the reference frame K^, 
rotating with the angular velocity Co, vectors e m and en are 
constant in time and the st ar rot a tes ar ound vector e m with 
the angular velocity (— fi p ) Link (2003), see fig. IAI1 where 

fi P = -e m fi tan(/3) (A6) 

Take, for example, any vector e a rotating together with the 
star as, see fig. I All 

Hp 

^ = [{1x4] (A7) 
at 

Hence in the reference frame K u it is possible to write 
de a - 

— — = [fix e a ] — [Co x e a ] = 
at 

= [Qp X 6 a ] 

Because of radio radiation may be generated at a some dis- 
tance from the center of pulsar tube, the precession may 
be manifeste d as subpulse s drift or variation of pulse profile 
lAsseol l|2002l ); iRudermanl l|200rj ). In order to exclude such 
manifestations it is necessary to assume that average distri- 
bution of radio sources in pulsar tube is axisymetrical. Also, 
it needs to assume that pulsar tube crossection is circular or, 
at least, its profile depend on vectors f2 and m and does not 
depend on small scale magnet ic field. The model of pulsar 
tube presented in [Biggs (1990) satisfies this criteria. 



